Electronic and magnetic properties of a hexanuclear ferric wheel 
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The electronic and magnetic properties of the hexanuclear ferric wheel [LiFee(OCH3)i2- 
(dbm)6]PF6 have been studied with all-electron Hartree-Fock and full-potential density functional 
calculations. The best agreement for the magnetic exchange coupling is at the level of the B3LYP hy- 
brid functional. Surprisingly, the Hartree-Fock approximation gives the wrong sign for the exchange 
coupling. The local density approximation and the gradient corrected functional PBE strongly over- 
estimate the exchange coupling due to the too large derealization of the d-orbitals. These findings 
are supported by results from the Mulliken population analysis for the magnetic moments and the 
charge on the individual atoms. 



I. INTRODUCTION 

Molecular magnetism has received great attention in 
the past few years, especially after a series of molecules 
has been discovered which show magnetism due to a 
molecular origin (see the review articles 0, Q)- One 
particular type of molecular magnets are the so-called 
ferric wheels where iron atoms are arranged ring-like. 
So far, wheels with 6 up to 18 iron atoms have been 
synthesized lj. This has also triggered interest among 
theorists. The description of magnetism still poses chal- 
lenges, and only in the last years ab-initio calculations of 
large molecular magnets have become feasible 0. 

The early studies employed wave- function based meth- 
ods: for example, the exchange interaction in transi- 
tion metal oxides was investigated at the level of peri- 
odic Hartree-Fock calculations Usually, it is found 
that the Hartree-Fock approximation strongly underesti- 
mates the exchange coupling (e.g. 0, [S E|) This was 
also demonstrated for smaller molecules, where demand- 
ing configuration interaction calculations were feasible: 
it was shown, that for a proper description of the mag- 
netic interactions, highly correlated wave functions are 
necessary 00- 

With the help of an appropriate embedding, the bulk 
could be approximatively described with clusters. Then, 
wave-function based correlation methods could be ap- 
plied and a value in reasonable agreement with the ex- 
perimental value of the exchange coupling was deduced 

%m 

In the past few years, density functional calculations 
have become more and more widespread, and have also 
been applied to molecular magnets (e.g. 0, 0, H^l. 
They have the advantage of including approximatively 
the effects of electronic correlation, at a very low com- 
putational expense. However, it is often observed 
that exchange couplings, computed at the density func- 
tional level (especially the local density approximation), 
strongly overestimate the experimental value. In a com- 
parative study, it was shown that hybrid functionals, 
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which include the exact (Fock) exchange, improve this 
and provide surprisingly accurate values 13] . This was 
also observed earlier in cluster calculations which were 
used to model the solid [tH ]. 

The ferric wheel with six iron atoms is one of the 
smaller molecular magnets and thus ab-initio calculations 
can be performed at a relatively low expense, besides 
the approach using model Hamiltonia ns [la, [la. Various 
molecules have been synthesized 0, Il8t Il9l I20L Ell I22I 
l23l l24l I25L I2H , [2tJ and techniques such as nuclear mag- 
netic resonance, high-field DC and pulsed-field differen- 
tial magnetization experiments, susceptibility measure- 
ments, high field torque magnetometry, inelastic neutron 
scattering and heat capacity measurements have been ap- 
plied. 

The experimental values for the exchange couplings, 
extracted from susceptibility measurements or from the 
sin glet -triplet gap, are typically around ^20 K [H I23I I24I 
I25L |27|. although also larger values have been suggested 
(e. g. 38 KI2CI). As there are still only few ab-initio 
calculations [28l l29| . these systems are therefore an inter- 
esting object to study. 

This article will deal with one of the six-membered 
iron ferric wheels [19| . In previous studies on a similar 
system with a gradient corrected functional, an overesti- 
mation of the calculated exchange coupling by a factor 
of ~ 4 was found and it was argued that schemes such as 
LDA+(7 or self-interaction correction might resolve this 
problem [2^, [2^|. As was argued above, hybrid function- 
als might be a different way to deal with this problem. 
We therefore studied this ring at the level of various func- 
tionals including hybrid functionals, and the traditional 
Hartree-Fock approach. 



II. METHOD 

All the calculations were done with the code CRYS- 
TAL2003 [ME!. We employed the Hartree-Fock (HF) 
method (in the variant of unrestricted Hartree-Fock), 
the local density approximation (LDA) with the corre- 
lation functional of Perdew and ZungerJ^, the gradi- 
ent corrected functional of Perdew, Burke and Ernzer- 
hof (PBE) [13 and the hybrid functional B3LYP (a func- 
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tional with admixtures, amongst others, of functionals by 
Becke, Lee, Yang and Parr). As described later on in this 
section, we performed the calculations mainly on a molec- 
ular system. Therefore, we also used a basis set library 
for molecular basis sets|34J rather then the CRYSTAL 
basis set library which contains essentially basis sets for 
periodic systems. We chose a [5s4p2d] basis set for iron 
|35| . where the diffuse exponent 0.041148 was omitted, 
a [2s] basis set for hydrogen^, a [3s2p] basis set for 
carbonj^ where the exponent 0.1687144 was replaced 
with 0.25 to avoid convergence problems, a [3s2p] basis 
set for oxvgen|38| with an additional d-exponent of 0.8 
and a [3s2p] basis set for lithium j^j where the exponent 
0.2 was omitted. After the mentioned modifications, the 
final basis sets were thus of the size [4s3p2d] (iron) , [5s^p] 
(carbon), [3s2pld] (oxygen), [2s] (hydrogen) and [2s lp] 
(lithium) . 

To explore the influence of further basis functions, vari- 
ous tests were performed: calculations with an additional 
sp-shell (exponent 0.04) and an additional d-shell (expo- 
nent 0.2) at the iron atom have been carried out, at the 
B3LYP and HF level. Also calculations with an addi- 
tional sp-shell (exponent 0.12) at the oxygen atoms were 
performed at the B3LYP level. This diffuse exponent led 
to linear dependence problems at the HF level. However, 
a HF calculation where the outermost exponent 0.270006 
for oxygen was replaced with two exponents (0.5, 0.2) was 
numerically stable. No significant changes for the com- 
puted exchange couplings were observed in any of these 
aforementioned tests. Therefore the basis sets, as de- 
scribed in the first paragraph of this section, can be con- 
sidered reliable, and all the results were obtained with 
these basis sets, if not explicitly stated otherwise. 

The properties of the ferromagnetic (FM) state (all 
spins parallel, total spin 30 [1b) and of the antiferromag- 
netic (AF) state (spins alternating up and down, total 
spin 0) were computed with each method. To obtain the 
local magnetic moments, a Mulliken population analysis 
was performed. 

The calculations were performed on an isolated 
molecule where the full symmetry of the molecule was 
exploited (ferromagnet: antiferromagnet: C3). The 
geometry was chosen according to the measurements by 
Abbati et al pjj and is displayed in figure [I] The iso- 
lated molecule was charged with +1, as lithium is ionized 
and charge neutrality is restored by a PFg ion for this 
particular ferric wheel. To assess the influence of the 
centered lithium ion it was also removed in some of the 
calculations. 



III. RESULTS 

The total energies, the difference in total energies of the 
ferromagnet and the antiferromagnet and the exchange 
parameters J are given in tabled To obtain the exchange 
parameters J from the calculation, the spin Hamiltonian 
of the Ising model was used: 
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H = J Si ■ Si+i — J ■ Se ■ Si 

i=l 

The variable i corresponds to the site index of the iron 
atoms, which in this model possess a spin of S = 5/2. 
With this Hamiltonian, the energy difference between 
ferromagnet (FM) and antiferromagnet (AF) then corre- 
sponds to: 

Ef ? - Eg = -12 JS 2 

Therefore, the exchange parameter J results from the 
difference of large numbers. Still the results are numeri- 
cally stable, as can be seen in table |U For the particular 
ferric wheel under consideration, the experimental value 
for the exchange coupling was determined to be -21 K 

am. 

This can now compared with the computed values. 
One interesting result is the ferromagnetic exchange pa- 
rameter obtained at the Hartree-Fock level (+7 K ) which 
does not reproduce the real antiferromagnetic nature of 
the molecule. This result was stable (±2K) with respect 
to all possible variations in the basis set as mentioned 
in section m and also when the central lithium ion was 
removed. 

Such a problem did not show up with any of the density 
functionals. The magnitude of J of the B3LYP calcula- 
tion (-31 K) is about 50% larger than the experimen- 
tal value. In the PBE and LDA calculations, it is even 
overestimated by a factor of ~ 5 (PBE: -108 K, LDA: 
-120 K). This agrees well with an earlier calculation for 
a similar ferric wheel by Postnikov et al [2^, [2(| with 
the PBE functional, where J was also largely overesti- 
mated (-80 K). The authors argued, that the d-orbitals 
in density functional calculations with functionals such 
as LDA or PBE, were not sufficiently localized to com- 
pare with experiment, and therefore the resulting ex- 
change parameters would overestimate the experimental 
values. They therefore suggested to use methods such 
as LDA+J7 or self-interaction corrections to improve this 
shortcoming. The B3LYP functional may be viewed as 
another possibility to remedy this problem, as it interpo- 
lates between the pure Hartree-Fock approach (with an 
unscreened Coulomb interaction and the correct Fock ex- 
change, but without electronic correlation), and standard 
functionals (including electronic correlation, but too de- 
localized d-orbitals). This was already demonstrated for 
bulk NiO where the B3LYP value for the exchange 
couplings was also found to be about 50 % larger than 
the experimental value, and LDA overestimated also by 
a factor of > 4. Also, the total density of states was 
calculated at the HF and B3LYP level and is shown in 
figure |2 Note that these calculations were carried out 
as bulk calculations. The HOMO-LUMO gap at the HF 
level is about ten times larger than at the B3LYP level, 
but there is no significant difference between the ferro- 
magnetic and the antiferromagnetic state at one level of 
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theory. However, the smaller gap will favor hybridization 
and thus the more delocalized situation. 

Comparing the results, the wrong sign for J at the HF 
level is rather surprising as it was usually observed that 
the Hartree-Fock result underestimated the experimental 
value, but with the correct sign, with few exceptions 
A possible source for the discrepancy can be found in 
the geometry of the molecule. For a system where the 
magnetic centers are in line with the bonding bridge 
atom, consideration of correlation leads to a strong in- 
crease of the antiferromagnetic exchangeparameter as it 
was found for NiO by de Graaf et al |!| and for other 
complexes by Fink et al0. In our case, the Fe-O-Fe 
bonding angle is in the range of 100 degrees^ and 
for such systems (and similarly for other ferric wheels, 
see e.g. Waldmann et al 39]) the exchange parameter is 
supposed to be less antiferromagnetic or even ferromag- 
netic compared to in-line geometries which is a result of 
the ordering of the magnetic orbitals according to the 
Goodenough-Kanamori rules [j^. In addition, due to the 
to large HOMO-LUMO gap, the system is too ionic at 
the HF level and thus the overlap of the orbitals is under- 
estimated, which also has an impact on the magnitude of 
the computed exchange interaction. Thus the small an- 
tiferromagnetic coupling and the underestimation of the 
exchange parameter due to the neglect of correlation at 
the HF level may be the reason for the wrong sign. 

The central lithium atom has an electrostatic influence 
to the molecule and the lithium ionic radius affects the di- 
ameter of the ferric wheel. Therefore, when the lithium is 
removed and the geometry input is not changed as in our 
case, no impact on the exchange parameter is expected. 

To analyze the magnetic states, the magnetic mo- 
ments were computed. They are displayed in table [H] 
A decrease of the local magnetic moment of the iron 
atoms from HF over B3LYP and PBE to LDA is observ- 
able, where the moment becomes more delocalized and 
is transferred to the surrounding oxygen atoms. This is 
particularly obvious in the ferromagnetic state. In the 
antiferromagnetic state, the iron magnetic moments are 
virtually the same as in the ferromagnetic state, apart 
from the sign. However, the bridge oxygens carry virtu- 
ally no magnetic moment, as there is an overlap of up 
and down spin density at these sites, originating from 
two neighboring irons, which as a whole cancels. For the 
same reason there is no magnetic moment at the centered 
lithium, and even in the ferromagnetic case the lithium 
spin is virtually zero. 

This feature is reflected in the spin-density plots in fig- 
ure [3] and in larger resolution in figure 21 For a more de- 
tailed study, an one-dimensional plot of the spin-density 
at the LDA level along a line from O (apical) via Fe, 
O (bridge) and Fe to O (apical) was performed (figure EJ). 
The iron atoms subsist in a high-spin state 3d^a^, but 
except for the HF method the moment is closer to 
as it was also obtained by Postnikov et al for a slightly 
different molecule [28l l29l |. As can be seen in table ITTI 
the local magnetic moment of the apical oxygen atoms 



increases from HF over B3LYP and PBE to LDA. In ad- 
dition, there is a local magnetic moment at the carbon 
C(l) atom (see figure ^ and table [HJ) which is about 4 
times (HF) to 2 times (LDA) smaller than the local mag- 
netic moment of the apical oxygen atom. Corresponding 
to the moment of the C(l) atom there is a small mag- 
netic moment at the carbon atoms C(2) and C(3), but 
with opposite sign and a much smaller value. Further on 
there is actually a moment at the first carbon atom (next 
to the carbon atoms C(2) and C(3)) of each CqH 5 ring. 
Its value is of the same order as the value of the C(2) and 
C(3) carbon atoms moments, but with opposite sign. In 
the ferromagnetic state, there is a small magnetic mo- 
ment at the C(4) bridge carbon atoms corresponding to 
the magnetic moment of the bridge oxygen atoms which 
is also, as already mentioned above, increasing from HF 
over B3LYP and PBE to LDA. 

The results from the Mulliken population analysis of 
the iron, oxygen (bridge and apical), lithium and carbon 
atoms are given in tables IIIII and IIVI The net charge 
of iron is between 1.27 and 2.16, at the various levels of 
theory, and thus far away from a formal charge 19J of +3. 
Again, Hartree-Fock gives a very localized picture with 
the largest net charge (2.16), and LDA a delocalized pic- 
ture with an iron charge of only 1.27. Oxygen carries 
essentially a single negative charge (~ -0.8 in LDA, ~ 
-1.1 in HF), the carbon atoms which have oxygen as a 
neighbor all have donated charge (0.4 - 0.8 |e|, depend- 
ing on the site and the level of theory). The lithium 
charge is between 0.4 (LDA) and 0.8 (HF). A plot of the 
difference of the charge density was performed between 
the values at the LDA and HF level (figure HJ|) . There 
is a positive charge density difference at the iron atoms 
because of the smaller charge at the LDA level compared 
to HF, the results for the other atoms follow from tables 
Iflll and IIVI This plot thus visualizes the more covalent 
picture obtained at the LDA level. The stronger cova- 
lency also becomes obvious in the Mulliken overlap pop- 
ulation: for example, the overlap population between Fe 
and the neighbouring oxygens atoms is 0.11-0.12 at the 
LDA level, but only 0.07-0.08 at the HF level. 



IV. SUMMARY 

Ab-initio calculations for a hexanuclear ferric wheel 
|l9| were performed. The exchange coupling parameter J 
was computed at various levels of theory, with the B3LYP 
functional providing the best agreement with experiment 
(-31 K, experiment: -21 K). LDA and PBE were found 
to grossly overestimate the exchange coupling parame- 
ter due to the too large delocalization of the d-orbitals. 
Surprisingly, the calculation at the Hartree-Fock level led 
to a ferromagnetic coupling. The electronic population 
and the spin densities were calculated. The total spin is 
distributed over several sites, besides the iron atom all of 
the neighboring oxygens carry some spin. The delocaliza- 
tion increases from HF over B3LYP, PBE to LDA which 
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shows the strongest derealization and thus the smallest 
magnetic moment on the iron site. 
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TABLE I: Total energies (Eh = Hartree), the differences in total energies for the ferromagnetic (FM) and the antiferromagnetic 
(AF) state and exchange parameters J 

method FM total energy (Eh) AF total energy (Eh) difference of total energy (mEh) exchange parameter J (K) 



HF -13295.7328681 -13295.7312480 -1.62 +7 

HF (without Li) -13287.8311111 -13287.8293232 -1.79 +8 

B3LYP -13334.5067552 -13334.5140876 7.33 -31 

PBE -13329.9510763 -13329.9768112 25.74 -108 

LDA -13273.9493819 -13273.9779682 28.58 -120 



TABLE II: Local magnetic moments at various sites, in /is, for the ferromagnetic (FM) and the antiferromagnetic (AF) state 
method Fe O (apical) O (bridge) C(l) C(2)/C(3) C(4) 





AF 


FM 


AF 


FM 


AF 


FM 


AF FM AF 


FM 


AF 


FM 


HF 


±4.79 


±4.79 


±0.04 


±0.04 


w 


±0.07 


±0.008 ±0.009 =f0.006 


-0.006 


w 


« 


B3LYP 


±4.32 


±4.34 


±0.10 


±0.10 


w 


±0.20 


±0.031 ±0.032 =f0.007 


-0.007 


w 


±0.005 


PBE 


±4.10 


±4.19 


±0.10 


±0.11 


w 


±0.22 


±0.041 ±0.046 =f0.005 


-0.003 


w 


±0.009 


LDA 


±3.96 


±3.94 


±0.11 


±0.11 


w 


±0.20 


±0.040 ±0.047 =f0.003 


-0.006 


w 


±0.010 



TABLE III: Mulliken charge of Fe, in \e\. Note that the charge is virtually identical for the ferromagnetic and the antiferro- 
magnetic state. 



method net charge s p d 



HF 
B3LYP 
PBE 
LDA 



2.16 
1.56 
1.40 
1.27 



6.28 12.28 5.29 

6.38 12.37 5.69 

6.39 12.37 5.84 
6.43 12.41 5.89 



TABLE IV: Charge at various sites, in |e|. Note that the charge is virtually identical for the ferromagnetic and the antiferro- 
magnetic state. 



method 


Li 


O (apical) 


O (bridge) 


C(l) 


C(2)/C(3) 


C(4) 


HF 


0.77 


-1.10 


-1.18 


-0.22 


0.82 


0.50 


B3LYP 


0.53 


-0.86 


-0.88 


-0.11 


0.62 


0.45 


PBE 


0.47 


-0.79 


-0.80 


-0.10 


0.56 


0.38 


LDA 


0.38 


-0.76 


-0.75 


-0.09 


0.56 


0.36 
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FIG. 2: The total density of states for the ferromagnetic and antiferromagnetic state at the B3LYP level (upper two graphs) 
and the HF level (lower two graphs). For this purpose, the calculations had to be carried out as bulk calculations, with a lattice 
constant of 100 A to avoid interactions between the molecules, and with the central lithium ion removed in order to have a 
non-charged unit cell. The top of the valence band is at -0.12 E h (B3LYP) and at -0.26 E h (HF). 
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FIG. 3: Spin densities at the LDA level for the ferromagnetic (FM) and the antiferromagnetic (AF) state. For all figures, the 
contour lines range from -0.0004 to 0.0005 in steps of 0.000035 electrons/(a.u.) 3 . Full lines indicate positive spin density and 
dashed lines indicate negative spin density. 
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FIG. 4: Rescaled plots of the spin densities for the ferromagnetic (FM) state at the LDA level (first graph) and for the 
antiferromagnetic (AF) state at the LDA (second graph) and the HF level (third graph) . 
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FIG. 5: One-dimensional plot of the spin density at the LDA level for the antiferromagnetic state along a line between certain 
atoms of the molecule, see labeling of the abscissa and insert caption. The scaling of the ordinate refers to the values in table 

un 




FIG. 6: Rescaled plot of the charge density difference between the LDA and the HF level for the antiferromagnetic state. The 
contour lines range from -0.005 to 0.005 in steps of 0.0007 electrons/ (a.u.) 3 . Full lines indicate positive charge density difference 
and dashed lines indicate negative charge density difference. 




